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Аналитическое решение контактной задачи о внедрении 
сферического индентора в мягкий упругий слой” 


С. С. Волков 
(Донской государственный технический университет) 


Приводится эффективный метод построения приближённого аналитического решения контактной задачи о 
внедрении сферического индентора в мягкий упругий слой, лежащий на жёстком упругом основании. Метод 
позволяет получить решение задачи, эффективное для всего диапазона значений характерных геометриче- 
ских и физических параметров задачи. Рассмотрены случаи отличия модулей Юнга слоя и подложки более, 
чем в 100 раз. Полученные контактные напряжения для различной степени жёсткости основания сравнива- 
ются с ранее известными результатами для случая недеформируемого основания; показано, что при отличии 
свойств слоя и подложки в 100 раз, максимальная разница между результатами составляет менее 3 %. 
Определена связь между вдавливающей силой и радиусом сферы. Задача поставлена в связи с проблемой 
контроля упругих свойств тонких упругих покрытий. 

Ключевые слова: контактная задача, сферический индентор, эффективное решение, мягкий упругий слой, 
упругое основание, модуль Юнга, контроль упругих свойств. 


Введение. Контактные задачи для упругого слоя на протяжении длительного периода времени 
привлекают внимание исследователей. Решения этих задач имеют многочисленные приложения, 
в частности, они эффективно используются при исследовании свойств износостойких и защитных 
покрытий. Одним из важных практических приложений является возможность интерпретировать 
результаты индентирования и наноиндентирования [1]. Результаты работы могут быть использо- 
ваны при изучении упругих свойств мягких материалов (мягкие полимеры, биологические ткани), 
когда следует учитывать влияние жёсткости подложки на образец. 

Решению контактных задач для упругого слоя посвящено большое количество работ, так, 
например, в [2] постановка классической задачи Герца обобщается, и рассматривается случай, 
когда сферический индентор вдавливается в слой, лежащий на недеформируемом основании или 
на однородном полупространстве (однородное полупространство имеет упругие свойства, отлич- 
ные от упругих свойств слоя). Задача сводится к решению интегрального уравнения. Приближён- 
ное решение интегрального уравнения строится двумя методами: сингулярным методом, который 
эффективен для очень тонких слоёв, лежащих на недеформируемом основании, когда толщина 
слоя Н меньше радиуса контакта а (^ =Н/а< 0:5) и регулярным методом, который применим в 


случае, когда толщина слоя такова, что ^>1 [3]. 

В данной работе для решения контактной задачи о внедрении сферического индентора в 
мягкий однородный слой на жёстком упругом полупространстве используется двусторонне- 
асимптотический метод [4, 5]. Под мягкостью слоя следует понимать то, что модуль Юнга рас- 
сматриваемого слоя в несколько раз меньше модуля Юнга подстилающего полупространства. За 
это отличие упругих модулей основания и слоя отвечает параметр относительной жёсткости ос- 
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нования п (относительной мягкости слоя), равный отношению соответствующих упругих модулей. 
Ниже рассмотрены случаи, когда параметр п = 2,5; 10; 100; 1000 и показано, что при увеличении 
значения этого параметра результаты, полученные двусторонне-асимптотическим методом, схо- 
дятся к известным решениям, полученным ранее [6] для недеформируемого основания. Исполь- 
зуемый метод позволяет получить эффективное решение задачи во всём диапазоне значений ха- 
рактерного геометрического параметра задачи Л. 





о АЕ 


Рис. 1. Постановка задачи 


Постановка задачи. Пусть сферический индентор вдавливается в поверхность мягкого упругого 
однородного слоя, сцепленного с упругим жёстким полупространством, силой Р (рис 1). С 
полупространством связана цилиндрическая система координат (Р7ф:Е Предполагается, что все 


деформации упруги и размер зоны контакта а мал по сравнению с радиусом К сферы, описы- 
вающей форму индентора. Силы трения между индентором и поверхностью полупространства 
отсутствуют. Считаем, что в окрестности начальной точки контакта форма сферического 
индентора аппроксимируется параболой 

=ф(г) =В/* (1) 
Аппроксимация обоснована для однородных тел [7] для малых радиусов контакта а < 0,1А . Это 
условие выполняется практически для всех случаев упругого сферического внедрения. Вне 
индентора поверхность полупространства не загружена. Под действием силы Р индентор 
перемещается на расстояние Хх вдоль оси 2. Коэффициенты Ламе изменяются с глубиной по 
следующему закону: 
1) ЛЕЛ,, М=М,, -Н<2<0; 
2) Л=Д, М=М,, -®<2<-Н. 
При сделанных предположениях граничные условия будут иметь вид: 


2=о,т. НО й 


и2=-Н 2 Оо = 
=-5(г)=-(х-4(г)), г <а о ии, и® и 


При (/; —2) > ® напряжения в полупространстве исчезают. Требуется определить перемещение 


(2) 


штампа и распределение контактных нормальных напряжений под штампом: 


90° |1 =-9 (г), г<а. (3) 
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Построение решения задачи. Используя метод интегральных преобразований, задача сводит- 
ся к решению следующего интегрального уравнения [4]: 


[®) (ура) о (и). (ирл" ди = @, (0) (г), г <1 (4) 


Здесь трансформанта (и) строится численно методом моделирующих функций [8]. 
В [4, 5] описан метод аппроксимации трансформанты ядра уравнения выражениями вида 
2 м 
(и) = (и) + (и) ‚где № (и яка ны (5) 
где А,В,,О, ЕС; С, ЕК — некоторые константы. 
Для трансформанты ядра (и) =/" в работе [4] получено аналитическое выражение, ко- 
торое служит для приближённого определения контактных напряжений и вдавливающей силы: 


<) =) ий ом сА» А «| ОРТ, (6) 
4аз© (0) 


ИЕ 1 (0+ >С, (®(А^)-А дтвен(4) | (7) 


Ю 1=1 


Неизвестные С, определяются из решения системы линейных алгебраических уравнений [4]. 


При рассмотрении контактной задачи о взаимодействии штампа с упругим слоем, лежа- 
щим на упругой подложке, в случае существенного отличия свойств слоя и подложки, использо- 
валось предположение, что слой лежит на недеформируемом основании [3, 6]. Для этой модели 
основания трансформанта ядра [(0) = 0 

Одним из условий использования двусторонне-асимптотического метода является то, что 
трансформанта ядра является всюду строго положительной, и более того, нигде не обращается в 
ноль. Предположение, что подложка под покрытием остаётся упругой, хотя её свойства отлича- 
ются от свойств покрытия в 100, 1000, и более раз, позволяет удовлетворить этому условию. 

Ниже считаем, что основание, на котором лежит слой, остаётся упругим, хотя его упругие 
свойства значительно отличаются от упругих свойств покрытия. 


(а) 
10. 
10 
10° 
10 
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Рис. 2. Трансформанты ядер интегрального уравнения. Кривые 1—5 соответствуют отличию модулей Юнга слоя и 
основания в 2,5; 10; 100; 1000 раз соответственно. Кривая 6 соответствует трансформанте ядра для 
недеформируемого основания 
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На рис. 2 изображены графики функций (и) — трансформант ядра интегрального урав- 


нения (4), построенные методом моделирующих функций [8]. Графики изображены в логарифми- 
ческой шкале по обеим осям. Кривые с номерами 1, 2, 3, 4, 5 соответствуют трансформантам ядра 
для материалов, у которых упругие свойства покрытия и подложки отличаются в 2, 5, 10, 100, 
1000 раз соответственно. Номером 6 обозначен график трансформанты ядра, для случая неде- 
формируемого основания, при построении которого была взята функция из [3], следующего вида: 
т 2х1 (2и) -4и 

й 2хсп(2и)+1+Х? +4и?' 
Анализируя рис. 2 можно увидеть, что по мере того, как увеличивается жёсткость основания, 
графики функций трансформант приближаются к своему предельному графику 6, для недефор- 
мируемого основания. Покажем, что значения контактных напряжений с увеличением жёсткости 
подложки приближаются к контактным напряжениям, соответствующим недеформируемой под- 
ложке. 
Численные примеры. Рассмотрим внедрение параболического индентора в однородный слой, 
лежащий на жёстком упругом основании. Параметр п, характеризующий относительную 


где х=3-4у. 





жёсткость, принимает значения 2,5; 10; 100; 1000 (п=Е, [Е (п))°). Коэффициент Пуассона, не 


нарушая общности, полагаем равным 0,3. 
Введём функцию Д, (и), характеризующую погрешность аппроксимации трансформан- 


ты ядра: 
А, (и) = (№ (и)- (и) (и) -100 %. (8) 
Формула (8) введена для оценки погрешности аппроксимации трансформанты ядра выражениями 
вида /\ (и). 
В таблице 1 приведены значения величины т“ = 2Ат (0) / ©, (0)а, характеризующей кон- 


тактные напряжения под штампом в точке г = 0. 



























































Таблица 1 
Контактные напряжения 
плах ^, (и), и = [0..5) 1,4 % 28% 3,3% 3% 3,2 % — 
№ 8 8 10 10 20 — 
п 2 5 10 100 1000 — 
Л = 0,25 4,098 6,375 7,991 10,075 10,326 10,342" 
лЛ=1 2,866 3,144 3,252 3,342 3,389 3,305' 
лЛ=4 2,562 2,555 2,560 2,556 2,542 2,568! 
В таблице 2 приведены значения величины Р” = 2КР/©, (0)а° ‚ определяющей связь меж- 

ду вдавливающей силой и радиусом сферы. 

Таблица 2 

Связь вдавливающей силы и осадки 

п 2 5 10 100 1000 — 
Л = 0,25 7,938 11,557 13,973 17,048 17,447 17,66' 
л=1 5,895 6,365 6,553 6,709 6,771 6,82' 
Лл=4 5,363 5,355 5,365 5,360 5,336 5,38! 























* Результаты из работы [6] для недеформируемого основания. 
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В таблице 1 и 2 значения крайнего правого столбца взяты из монографии [6] и характери- 
зуют т’и Р” для случая вдавливания параболического штампа в упругий слой, лежащий на не- 
деформируемом основании. 

Из анализа таблицы 1 видно, что максимальная разница между контактными напряжения- 

ми при п = 1000 и значениями для недеформируемого основания не превышает 3 % при Л = 1. 
Заметим, что из таблиц 1, 2 видно, что значения т’и Р” принимают близкие значения при 
П = 100 ип = 1000. Обратим внимание также на то, что упругие свойства основания не оказыва- 
ют сильного влияния на распределение контактных напряжений, для значений параметра Л > 1 
(если сравнивать значения напряжений при п, равном 2 и 1000, то максимальная разница не пре- 
вышает 15 % для Л = 1). 
Заключение. Использование двусторонне-асимптотического метода при решении осесиммет- 
ричной контактной задачи о мягком упругом слое на жёстком основании позволяет учесть дефор- 
мируемость подложки под мягким покрытием твёрдого основания и получить решение, эффек- 
тивное для всего диапазона значений характерных геометрических и физических параметров за- 
дачи, что является важным при исследовании механических свойств биологических тканей и мяг- 
ких полимеров. 

Автор благодарит С. М. Айзиковича за постановку задачи и советы при реализации 
решения. 

Библиографический список 

1. Тэи!, У. Т. бирзгае еНес5 оп папотаепаНоп тесНапса! ргорейу теазигетепе оЁ ой 
Нип$ оп Вага зиб$га(ес / Т. У. Ти, С. М. РВагг // Зоигпа! оЕ Маепа!5 КезеагсН. — 1999. — Мо. 14. — 
№. 1. — Рр. 292—301. 

2. Е-ЗЛегЫтеу, М. С. 0. ТВе Негдап Сог{асЕ ог биасез Сомегед мИй МеаШс ЕИтз$ / 
М. С. О. Е-бВегытеу, 2. Напа // Меаг. — 1996. — Ма. 40. — №. 3. — Рр. 325—337. 

3. Ворович, И. И. Неклассические смешанные задачи теории упругости / И. И. Ворович, 
В. А. Александров, В. А. Бабешко. — Москва: Наука, 1974. — 456 с. 

4. АлкомсИ, 5. Апа!УЧса! зоиЧоп ог {Не эрпейса! паегпаНоп ргоЫет Гог а На[Н-зрасе ми 
дгаЧег{5 ми Че дер ва$Ис ргорег@ез / 5. Амомси, \. Аехапаго\, 3. КаКег, |. Кгепеу // Ттегпа- 
Чопа! Зоигпа! ог 5014$ апа 5гисигез 39. — 2002. — Рр. 2745—2772. 

5. АШКОМСИ, 5. ЕмааНоп ог {Пе вас ргорег4ез оЁ а РипсНопа!у-дгадеЯ соаЧта Нот Не т- 
ег(аНоп теазигетепё / $. АКомсН, |. Кгепем, 1. бемо$Напом, Т. Тибо // ХАММ. — 2011. — 
Рр. 1—23. 

6. Александров, В. М. Неклассические пространственные задачи механики контактных 
взаимодействий упругих тел / В. М. Александров, Д. А. Пожарский. — Москва: Факториал, 
1998. — 288 с. 

7. бигебй, 5. А пем тео Гог езИта Ипа гезача! $гез5ез Бу пугитещеа эВагр таепайоп / 
5. бигезй, А. Е. СаппаКороцо$ // Аа таег. — \а1. 46. — №. 16. — 1998. — Рр. 5755—5767. 

8. Бабешко, В. А. Методы построения матриц Грина для стратифицированного упругого 
полупространства / В. А. Бабешко, Е. В. Глушков, Н. В. Глушкова // Журнал вычислительной 
математики и математической физики. — 1987. — Т. 27. — № 1. — С. 93—101. 


Материал поступил в редакцию 09.08.2012. 
Вегегепсе$ 


1. Ти, У.Т., РВаг; С.М. зибугае еКес5 оп папотаеп(аНоп теспаптса! ргорецу теазигетепе 
ОЕ 50 ЕИт5 оп Вага зиб$гаез. Лоигпа/ о! Маепа Кежеагсй, 1999, \о|. 14, по. 1, рр. 292-301. 


Физико-математические науки 








2. Е-бНегытеу, М.С.0., Напд, 2. ТВе Нейжап Согкасе оЁ биПасез Соуегед МИ Меа!с ЕИтз. 
И/еаг 1996, \о1. 40, по. 3, рр. 325-337. 

3. Могомсй, Т.Т., Аехапагом, М.А., ВаБезИКо, М.А. Меказ/срезКуе зтезваппууе гайас!! ео! 
иргидо$й. [Мопбаз$!са! птхей ргоШеглз ог е!а$Нску {Пеогу.] Мозсоми: МаикКа, 1974, 456 р. (т Кизчап). 

4. АкомсИ, 5., Аехапагом, \., КаКег, )., Кгепем, |. Апа!уЧса! зо!иНоп о {Пе эрпейса! паега- 
Чоп ргоМет Юг а Па№-зрасе мА дгаФеп{ мИИ Пе дер еа$@с ргорегез. ГуЁетайопа! Лоигпа! оЁ 
50/45 апа ЗгисЁиге$, 2002, 39, рр. 2745-2772. 

5. АКОМСИ, 5., Кгепем, |., ЗемозНапом, Т., ТгибсМК, Т. Еуа!цаНоп ог {Ле ва$Ис ргореез оЁ а 
ГипсбопаЙу-дгаде4 соайпа гот {Ие таета{Ноп теазигетет5. ХАММ, 2011, рр. 1-23. 

6. Аехапагом, \.М., РогПаг$Юу, О.А. Мек/аз/све$Куе ргозвгап{уеппууе гадас!й! текваги Коп- 
{аКтпукВ угайтодеу$/у иргидИКВ {е/. [Мопааззса! зраНа! ргоМепз оЁ теспапт!с$ оЁ е!азЧс Бофез соп- 
{асЕ ищегасЧоп.] Мозсоми: РаКопа!, 1998, 288 р. (т Визчап). 

7. ЗигезН, 5., баппаКороЧ°$, А.Е. А пем теёо4 Юг езЧта{птд гез@чца! $гез5ез Бу тги- 
тепёе4 зПагр таегтаНоп. Аса таг 1998, уч. 46, по. 16, рр. 5755-5767. 

8. ВаБезНКо, М.А., Си$НКом, Е.М, Си$НКоуа, М.М. Меюду розбоуетуа та Сппа Ауа ${гай- 
Пёгоуапподо иргидодо ро/иргогап$уа. [Сбгееп тасез дечдп теоб$ Гог е!а5Чс угаНЙаЫе Па!- 
зрасе.] Хригпа/ иус/5Шепоу таетабй К! / таетайсйезКоу ИА, 1987, мо|. 27, по. 1, рр. 93-101 (т 
Кизчап). 


АМАГУТТСАЕ 5ОГОТТОМ ТО СОМТАСТ РКОВЕЕМ ОМ $РНЕКТСА! ТМОЕМТЕК 
РЕМЕТКАТТОМ Т1МТО $ОЕТ ЕГАЗТТС ГАУЕВ” 


$. $. МоКо\ 
(Роп За Тесптка! Упмег$Ку) 


Ап ейс/епЕ теоа Гог соп$ёгисипа ап арргохтаЁе апа/У&са! зо/ивоп 0 {пе сотасЕ ргоМет о! те репеёаНоп оЕ 
а зрвейса! паепЕег тю {Пе 5ой. е/а$8с !ауег утд оп а поа е/а$8с ЮипааНоп [5 ргорозеа. Те тепо4 реттй ю 
оБат йе роет зо/ийоп ейесвуе Гог те Пи! гапде о! пе ргоМет свагасейНс деотеё"са! апа рпуза! ра- 
гатеег иа/иез. боте сазез ипеп Уоипд’5 тоди/ о! {ве !ауЕг апа {те зибзгае (ег тоге тап а Випагедо!А аге 
апа/угеа. Тре офатеа сопасЕ з{теззез Юг ИТегепЕ уаиез о! {те ЮипааНоп 5$#Йте55 аге сотрагеа ийй {те Ю!- 
тепу Кпоит гези!5 Гог пе попаеюгта Ме Юипаайоп сазе. 1+ [5 пои паЁ ивеп #1е еазйс ргорегйех ог! те |ву- 
ег апа те зибзгайе ег а ВипагедаюЮ/, тахтит @#егепсе Бенмееп {теже гезий [5 [е55 пап 3%. Тре ге/айоп- 
5/р Бевиееп 1е ргезхтд Гогсе апа пе габ!из оЁ {Те зрПеге 1/5 а егттеа. Те диезйВоп 15 риЁ оилпд ю {те соп- 
о! ргоМет ог пе т соайпд$ е/а$8с ргорегИез. 

Кеуигога5: сошасЕ ргоМет, зррепса/ паетег ейесйуе зо/ибоп, 50й еа5йс ([вуер еазЁс Гоипдавоп, Уоипа5 
тоди/из, сопго! оЁ е/азНс ргорегйес, 


* Тре гееагсн {5 допе МИН Ше рагЧа! Япапса! зирро" гот ВЕЕТ (11-08-91168-СЕЕМ_а, 12-07-00639_а) апа от ТНе Митену 
ог ЕдисаНоп апа 5аепсе оЁ {Пе Визчап РедегаНоп (5{ае СопёгасЕ №. 11.519.11.3028, Адгеетег5 №5. 14.В37.21.1131, 
14.132.21.1693). 
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